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I. INTRODUCTION 

Short and intense laser pulses allow nowadays to drive 
atoms and molecules in nonperturbative regimes going 
from adiabatic (nano- and picosecond) to sudden or im- 
pulsive (femtosecond). Recent examples concern the 
alignment of molecules, which can be achieved during 
nanosecond pulses or after femtosecond pulses pj. Cor- 
rections to perfect adiabaticity can be analyzed in terms 
of superadiabaticjp, 13,0,01 and Davis-Dykhne-Pechukas 
techniques 0, Q, |J. On the opposite side, regimes be- 
yond the impulsive approximation, i.e. beyond the limit- 
ing case of pulses described as S— kicks, have not been yet 
much explored due to a lack of adapted tools of analysis. 

It is well-known that one can treat periodic perturba- 
tions using extended Hilbert spaces where time is con- 
sidered as a new dynamical variable, in order to render 
the problem autonomous |?j.llO(. This app roach, that can 
be formulated as Floquet theory 0,EaEl, allows one 
to eliminate systematically secular terms (i.e. terms that 
grow arbitrary with time), which would otherwise lead to 
divergences. Pulse-driven dynamics associated to Hamil- 
tonians localized in time requires a different treatment of 
secular terms. In this case, since the perturbation acts 
only during a finite time interval, the secular terms do 
not lead to divergences. 

This article contributes to develop a time-dependent 
perturbation technique, that is in particular suited for 
pulse-driven dynamics, on the basis of Refs. 0, . In 
Ref. [T^j. we constructed a superexponential perturba- 
tion theory which preserves the unitarity of the evolu- 
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tion operator at each order, and applied it beyond the 
impulsive regime by considering an expansion where the 
perturbative parameter is the characteristic duration of 
the time-dependent interaction compared to the charac- 
teristic time for the free evolution. We have shown that 
it converges in any regime (from impulsive to adiabatic) 
in two-level systems. This derivation is based on the 
Kolmogorov-Arnold-Moser (KA M) tec hniq ue applied in 
an extended Hilbert space [II El El El 13 • In Ref - 
|15| , we presented an improvment of this technique tak- 
ing advantage of free parameters, connected to secular 
terms, that are available to reduce the error without prior 
knowledge of the exact solution. This optimization en- 
hances the accuracy of the method in such a way that the 
first order approximation gives a satisfactory description 
up to fairly large values of the perturbative parameter. 

The present article contains a detailed description of 
the methods announced in |15|. Instead of using an ex- 
tended space, we formulate the derivation in a simpler 
way, by stating the perturbation iterations directly at 
the level of the evolution operator in the original Hilbert 
space. This scheme allows us to consider and compare in 
a unified way various time-dependent perturbation tech- 
niques. In particular we make the connection with the 
well-known Magnus expansion |2l|, that has been used 
by Henriksen et al. to construct an improved impulsive 
approximation (2^]. We also develop and investigate the 
accuracy optimization which can be applied to the time- 
dependent Poincare-Von Zeipel, the time-dependent Van 
Vleck and the time-dependent KAM techniques. 

The paper is organized as follows. In Sec. [H] we recall 
the Magnus expansion and outline the time-dependent 
versions of the Poincare-Von Zeipel, the Van Vleck and 
the KAM techniques. We highlight the free parameters 
and free operators that may be present in these unitary 
perturbative methods. In Sec. IIIII we exploit these de- 
grees of freedom to improve the accuracy of a given level 
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of approximation. Section IIVI is devoted to the applica- 
tion of these techniques beyond the impulsive regime and 
the illustration on a pulse-driven two-level system. The 
conclusions are given in Sec. and some details of the 
calculations are reported in Appendixes A to C. 

II. UNITARY TIME-DEPENDENT 
PERTURBATION THEORIES 

We consider the Schrodinger equation 

i%-U Hl (t,to)=H 1 (t)UH 1 (t,t ), U H dto,t ) = t, (la) 
at 

where H\(t) is a time-dependent matrix or a time- 
dependent operator in a Hilbert space H. Assuming that 
one can decompose H\ (t) according to 

H 1 (t)=H (t)+eV 1 (t), (lb) 

where H (t) is such that its propagator E/y (i,io) is 
known and V\(t) is localized in time (i.e., vanishes outside 
a finite interval), we are looking for a unitary perturba- 
tive expansion of the full propagator Ujji (t, to)- This can 
be achieved by two classes of techniques that we outline 
below: i) the order by order methods, namely the Mag- 
nus expansion, the time-dependent Poincare-Von Zeipel 
technique and the time-dependent Van Vleck technique, 
where after n steps the remainder is of order e n+ , and 
ii) the surperexponential KAM technique where the re- 
mainder is of order e 2 

Below we will have to consider the propagator 
Uuxitito) in the interaction representation with respect 
to Ho(t) 

U Hii (t,t ;s) = UH (s,t)U Hl (t,t )U Ho (to,s), (2a) 



where s is an arbitrary time (the standard interaction 
representation corresponds to the case s — to). This 
propagator satisfies the Schrodinger equation 

i^ t U H i(t,to]s)=Hi(t,s)U Hi {t,t -,s), (2b) 

and the associated Hamiltonian reads 

Hi(t;s) = eU Ho (s,t)V 1 (t)U Ho (t,s). (2c) 

We will also consider a new representation defined with 
the help of a unitary transformation T(t; s) according to 

U„(t, t ; s) = Tt(t; s)U Hl (t, to)T(to, s), (3) 

where H(t;s) is a new Hamiltonian. This expression is 
reminiscent of Eq. I|2a|) although T(t; s) need not be a 
propagator but a unitary operator which features an ar- 
bitrary parameter s and satisfies the property T'(t\ s) = 
T(s;t). 

A. Magnus expansion 

The solution to Eq. (|la|) can always be put in the form 
of an exponential 

U Hl (t,t ) = e -^"i (**>), M Hl {t ;to) = 0, (4) 

where M^i^to) is self-adjoint to ensure unitarity. The 
exponent Mn^t^to) is generally not simply the integral 
of Hi(t) owing to the non-commutativity of this latter 
for different times. Indeed from Eq. Ijlaf) one deduces the 
following equation for Mjj 1 (t; to) 



j t M Hl (t;to) = + l - [M Hl (t; t ), #!(£)] - i [M Hl (t;t ), [M Hl (t;t ), H^t)}} + 



(5) 



We refer to the paper of Magnus |2lJ for a derivation of this equation (see also Ref. |2ij). The matrix or the operator 
Mfli(t; ^o) is obtained by integrating the first term on the right hand side of Eq. JSJ and substituting the result into 
the next terms of this equation. One then repeats this procedure, known as Picard's iteration, with the resulting 
second and subsequent terms. This yields the Magnus expansion 



M Hl (t;t„) 



/"* i /•* 

/ duH\(u) + - I du 
J to 2 J t(j 



dvH-i(v),Hi(u) 



du 



dv 



dwH^w),^^) 



dvHi(v), 



dwHi(w), Hi(u) 
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The number of terms in this expansion grows very rapidly 
[24|. Notice that this expansion is not limited to per- 
turbation theory [ i.e. to a Hamiltonian of the form of 



Eq. Ijlbfl ] although it is of interest only when the subse- 
quent terms are negligible. This will be the case if one is 
interested in obtaining an expression valid for very short 
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times only or in the presence of another small parameter. 

In the framework of perturbation theory, if we were to 
substitute Eq. (|lb(l into Eq. JJjJ we would generally obtain 
contributions to a given order in e from an infinite series 
of terms of the Magnus expansion |3^| . This can easily 
be circumvented by considering Eq. J5J, the equivalent 
of Eq. Q in the interaction representation. The Magnus 
expansion pertaining to Eq. (|2b|l allows one to write its 
solution in the form 

U H i (t, to; s) = exp j -ieM H{ (t; to, s) \ , (7) 

with eMyi (t; to, s) given by Eq. © where Hi(t) is re- 
placed by H\{u; s). By virtue of Eq. l|2c|l . each H\(u; s) 
carries a prefactor e. Hence, there is no e-independcnt 
term and terms of M H ^(t;t a ,s) with higher numbers of 

H\ (u; s) are now of higher orders in e. In order to display 
the explicit dependence on the parameter s, we decom- 
pose the propagator entering through H\(u; s) according 
to Uh (s,u) = UH a {s,t)UH {t,u) and note that the left- 
most UH (s,t) and rightmost UH (t,s) of each term of 
M ff i (t; to, s) factor out (while the inner ones cancel each 
other) 

oo 

eM H{ {t;to,s) = U„ (s,t)J2^M k {t;to)U„ {t,s). (8) 

fe=i 

Returning to the original representation with the help of 
Eq. ijSaJ). Eq. lj7|) yields 

U Hl (t,t ) = exp j^^e fc M fc (t;t )j U Ho (t,t ), (9) 

where use was made of the identity 

Ae B A~ 1 =e ABA ~\ (10) 

By truncating the infinite series of Eq. @ to order n 
one obtains the unitary approximation 

U Hl (t,t ) = U$(t,t ) + O(e n+1 ), (11) 
with the Magnus expansion 

Ujg (t, t a ) = exp | -i ^ M k(tl *o) | U Ho (t, to). (12) 

We stress that this expression is independent of the ar- 
bitrary time s chosen in Eq. [34| ■ Each term of the 
exponential can be cast into the form 

M k (t;t ) = f duU Ha (t,u)V k (u)U Ho (u,t), (13) 

where V k (u) is deduced from Eqs. (|2"c|l . (JHJ and (JSJ). The 

case k = 1 features the perturbation itself while for the 
first few values of k that we shall use below, dropping the 
time arguments, one finds 

V 2 = - l -[M 1 ,V 1 ], (14a) 
V s = -i [M 2 , Vi] + ~ [Mi, [Mr, Vi]] . (14b) 



B. Time-dependent Poincare-Von Zeipel expansion 

The classical Poincare-Von Zeipel technique was 
adapted to quantum mechanics by Scherer to treat both 
time-independent [2{| and time-dependent Ispft sytems. 
In the time-independent case it was shown [25j that this 
technique coincides with the usual Rayleigh-Schrodinger 
expansion. We present in Appendix 1X1 a time-dependent 
Poincare-Von Zeipel technique that has the advantages 
to be strictly unitary upon truncation at any given order 
and does not require the consideration of a so-called ex- 
tended Hilbert space. This method amounts to map the 
full propagator Uh^ it, to) into a new effective propagator 
Uh b (t, to) with the help of a single unitary transformation 
according to an equation similar to Eq. ©. Moreover, 
this formulation exhibits free parameters available to im- 
prove the accuracy of a given step of the algorithm by 
a procedure we describe in Sec. IIII Bl Finally, we show 
that the time-dependent Poincare-Von Zeipel method in- 
cludes the Magnus expansion as a particular case. 

C. Time-dependent Van Vleck technique 

To provide a general perspective of unitary time- 
dependent perturbation theories we introduce in Ap- 
pendix a time-dependent version of the Van Vleck 
technique that is more widely used than the preceding 
method in the stationary case |3l| . It consists in trans- 
forming the full propagator into a new effective propaga- 
tor iteratively through a series of unitary operations. Our 
purpose is actually not to introduce yet another variant 
of perturbation theory but to emphasize that this time- 
dependent version of a well-known technique is i) compa- 
rable to the Poincare-Von Zeipel method which we show 
to be closely related to the Magnus expansion, and ii) 
not as performant as the KAM technique detailed below 
since it is an order by order perturbative method. 

D. Time-dependent Kolmogorov-Arnold-Moser 
expansion 

The time-dependent KAM technique aims at obtain- 
ing a superexponcntial expansion for the propagator 
Unxitjto) through a series of unitary transformations. 
It is sometimes said to be superconvergent. However, 
superexponcntial is more appropriate since the depen- 
dence on e of the remainder after n iterations is indeed 
the exponential of an exponential (e 2 ) while the ac- 
tual convergence of the algorithm has to be examined 
specifically. 

1. First iteration 

The first step is to construct a unitary operator T\(t) 
which transforms the propagator Uh ± {t, to) we are look- 
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ing for into the propagator UH 2 (t,to) [cf. Eq. (J3J)] 

Tl{t)U Hl {t,t )Ti(t Q ) = U H2 (t,to), (15) 

where Ur 2 *o) is associated with the sum of an effective 
Hamiltonian Hf(t) which contains all contributions up to 
order e and a remainder e 2 V 2 (t) 



H 2 (t) = Hf(t)+e 2 V 2 (t). 



(16) 



As the new propagator is generated by a sum of Hamil- 
tonians it can be expressed in terms of the propagator 
Uii«(t, to) and a propagator Up 2 (t, to) by considering the 
interaction representation with respect to Hf(t;ti) 

Uh 2 (t, t ) = Uh} (t, t'{)U R2 (t, to; t'l)U Hl (t'l, to). (17) 

By virtue of Eq. the propagator Up 2 (t, to; t") satisfies 
a Schrodinger equation whose Hamiltonian is 



e 2 R 2 (t;t'{) = e 2 U H? (t'{,t)V 2 (t)U H c(t,t'{) 



(18) 



Being associated to a Hamiltonian of second order in e, 
this propagator will be neglected in Eq. (|T7Jl . i.e., re- 
placed by the identity. 

The effective Hamiltonian can be decomposed as 

flJ(t) = flo(t) + cI>i(t), (19) 
which allows one to express Unfit, to) as 

U H} (t,t ) = EW*,*i)Cft(*.*o5*i)£W*i,*o), (20) 

where Up^t^to) is the propagator corresponding to the 
Hamiltonian 



(21) 



Thus far the only restriction on the Hamiltonian tD\{t) 
is that it be of order e. Hence we have the freedom to 
choose it so as to be able to determine explicitly the 
propagator Up x (t,to;ti). This will be the case if 

Dx(t) = U Ha (t, ii)-Di(ti;ti)£/ Ho (ti,t) = Dx{t; h), (22) 

with Di{t\;t\) arbitrary. In this paper, we consider ex- 
plicitly two possibilities 



DxftijtiJsO, 



(23a) 
(23b) 



The first one is a trivial choice which gives nevertheless a 
nontrivial one- iteration KAM expansion (it will be shown 
in Sec. Ill D"4l to coincide, for the first iteration only, with 
the first order Magnus expansion). The choice to relate 
Di{ti;t\) to the perturbation according to Eq. (|23b|) is 
also rather natural as this operator enters the effective 
Hamiltonian (we shall discuss the fact that the pertur- 
bation is evaluated at an arbitrary time t\ in Sec. lIII A*t . 
Choosing one of the possibilities of Eq. (|23|1 implies that 
the propagator defined in Eq. (|2U|I reads 



The effective propagator can then be given a convenient 
form using Eqs. J20J) and JUJ 

U„ t (t,t Q ) = [/ Ho (Mo)e- l(t - to)£Dl( * o;tl) 

= U H ^t,t ;tx). (25) 

We now express the requirement that T\{t) defined in 
Eq. (|15fl be such that e 2 V 2 {t) contain no terms of order 
lower than e 2 . We first multiply Eq. (|lafl from the left by 
T±(t) and from the right by Ti(to) to deduce employing 
also Eq. l(T5jl 

e 2 V 2 (t) = Tlm^T^t) Ht(t) T 1 t (^)^|T 1 (^). (26) 
Writing T\ (t) in the exponential form 

Ti(i) = e- fcWi W, (27) 

we then require that all terms of order lower than e 2 in 
Eq. H26|) vanish identically. This leads to a differential 
equation for the self-adjoint operator W\(t) 

^W 1 (t)=V 1 (t)-D 1 (t;t 1 )+i[W 1 (t),H (t)}, (28) 

and defines the remainder e 2 V 2 (t) as the right side of 
Eq. |(2SJl. We stress that to arrive at Eq. we do not 
identify terms order by order. Hence this equation is still 
valid if any of the operator featured contains a further 
dependence on e. The general solution to Eq. I|28(l reads 

Wx{t)= f duU Ho (t,u)\yi(u)-Di(u-,ti)]UH (u,t) 

+ U Ho {t 1 to)B 1 U HQ {to 1 t), (29) 

where B± is any constant self-adjoint operator. In the 
present work we choose B\ = 0. 

Substituting Eq. fTTjl into Eq. (fT5|l and replacing 
UR 2 (t,to;t'{) by the identity as discussed above, one ob- 
tains the KAM approximation 



U Hl (t,t ) = U i H ) 1 (t,to) + 0(e 2 ), 
with the one-iteration KAM expansion 

Ugfato) = T^UH^to^lito). 



(30) 



(31) 



Up 1 (t,to;ti 



-i(t-to)eDi(ti;ti) 



(24) 



We emphasize the dependence of the following operators 
on the arbitrary times t±, t[ and t'{ 

H t{t) = Hf(t;ti), 
TiW^Ttf; h,t[), 
Wx{t) = W^tiJj), 
V^-V^t-M,^), 
R 2 (t;t'{) = R 2 (t;t u t[,t'{). (32) 

As a consequence, the one-iteration KAM expansion de- 
pends on t\ and t[ . In Scc. lIII Bl we shall show how these 
parameters can be chosen to improve the accuracy of the 
algorithm. In addition, recall that there are two con- 
stant operators, in Eq. J22J an d B 1 in Eq. |2HJl, 
that can be freely chosen. Note that with the choice of 
Eq. I|23a|) there is no dependence on t\. 
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2. Second iteration 

In the first iteration of the KAM algorithm we started 
with the Hamiltonian H\(t) = H (t) + eVi(t) and the 
known propagator UH (t,t ). We constructed, with the 
help ofTi(t), anew Hamiltonian H 2 (f) = Hf(t) + e 2 V 2 (t) 
and its propagator Un 2 (t, to). We approximated these 
operators to first order by retaining only the effective 
Hamiltonian Hf(t) and its propagator i ), dis- 

carding thus the remainder e 2 V2(t) and the related prop- 
agator U R2 (t,t ;t'{), 

To go one step further, unlike standard perturbation 
theory which reduces the size of the remainder from e 2 
to e 3 , the KAM algorithm takes advantage of the fact 
that after one iteration it produces a new perturbation 
e 2 V2(t) whose order is the square of that of the origi- 
nal perturbation eV\ (t) . Hence by considering H2 (t) and 
in particular the perturbation e 2 V2(t) as the new start- 
ing point, a KAM transformation T2(t) produces a new 
perturbation e^V^t) [whose order is indeed (e 2 ) 2 ]- This 
is the essence of the superexponential character of the 
KAM algorithm. The fact that the new perturbation is 
of order e 4 instead of e 3 as in a standard perturbation 
theory, allows one to anticipate the importance of keep- 
ing higher order terms in e 2 V2(t), in particular terms of 
order e 3 which would otherwise be absent if one further 
iterates the algorithm. 

The second KAM iteration amounts thus to reproduce 
the first iteration on the newly constructed Hamiltonian 
H 2 (t) = Hl(t) + e 2 V 2 {t): the effective Hamiltonian Hf(t) 
and propagator Unfit, to) now play the role of the pre- 
vious unperturbed Hamiltonian and propagator respec- 
tively, we replace e by e 2 and increase each subscript by 
one unit. The perturbation e 2 V2(t) is given by Eq. 
which we expand using Eqs. H27(l - (|28(l and the Hausdorff 
formula HtH 



1 



1 



e A Be- A = B + -[AB] + -[A 1 [A 1 B]] 



2! 



3! 



[A, [A, [A, B]]] 



(33) 



to obtain 



e 2 V 2 



Vi + D{\ 



-[^,[1^,2^ 




Di]} 
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[W x , [W 1 , [Wi,3Vi + Di] 



(34) 



We can rewrite this expression in a compact form with 
the shorthand notation ad fc (A, B) for the k nested com- 
mutators [A, ■ ■ ■ [A, [A, B}]-- ■} 



k = 
k > 1. 



(35) 



The perturbation we shall start from in the second iter- 
ation of the KAM algorithm reads thus 



e 2 V 2 



°° ,-fcJc+l 

E 1 ^ 



ad k (W 1 ,kV 1 +D 1 ). (36) 



We recall that the presence of a power series in e con- 
stitutes no difficulty for this algorithm. It is on the con- 
trary crucial in this superexponential technique to keep 
terms up to the final order one is interested in. Terms 
of order e 3 , for instance, do not appear through e 4 V3 (t) 
but through the second term in the series of Eq. (|34H or 
Eq. (U for e 2 V 2 (t). 

We now proceed to the construction of the second 
KAM iteration as described above. The unitary trans- 
formation T2(t) is such that 



T^t)U H2 (t,to)T 2 {to) = U H3 (t,t ), 



(37) 



where the propagator UH 3 {t,to) is associated with the 
Hamiltonian 



H 3 (t)^H c 2 (t) + e 4 V 3 (t), 
and can be expressed as 



(38) 



U H3 (t,t ) - U H5 (t,4)U Rn (t,h]1%)UH*(&,to). (39) 

The new effective Hamiltonian can be decomposed as 

H c 2 (t)^H c 1 (t)+e 2 D 2 (t), (40) 

and its propagator accordingly written in the form 

U H -{t,t ) = U m (t,t2)Up 2 (t,to;t2)U Hi (t2,to)- (41) 

From Eq. (j2J one deduces that Up 2 (t,to',t 2 ) is associated 
with the Hamiltonian 

e 2 P 2 (t;t 2 ) = U Hl {t2,t)D 2 {t)U H «{t,t 2 ), (42) 

where Unfit, to) is given by Eq. I|25(l . The correspond- 
ing Schrodinger equation is straightforwardly integrated 
if D2(t) is taken as 

D 2 (t) = U H i(t,t2;ti)D 2 (t 2 ;t 2 )U H} (t2,t;ti) 

= D 2 (t;t 2 ), (43) 



with D2 (t2 ; t2 ) arbitrary. Two appealing cases are 

D 2 (t 2 ;t 2 )=0, 
D 2 {t2;t 2 ) = V 2 (t 2 ). 

One then obtains 

Up 2 (t,to;t 2 ) ^ e-^' t0 ^ 2D2(t2 ' t2 l 



(44a) 
(44b) 



(45) 



The new effective propagator can be rewritten using 
Eq. m as 

U H «(t,t ) = J7 Ho (t^ ) e - l ( t -*«) e - D i( t O' t i)e- J (*" tt ») £2 - D2 ( to; * 2 ) 
= U H ^t,to;t 1 ,t 2 ). (46) 

We now come to the definition of T 2 (t) given in Eq. I|37|l 
and require that e A V 3 {t) contain no term of order lower 
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than e 4 . Inserting T 2 (t) into the Schrodinger equation for 
UH 3 (t,to) one arrives at 

eV 3 (t) = TKt)H 2 {t)T 2 {t) fff (t) Tl{t)i^T 2 {t). (47) 



We write 



T 2 (t) = e- ie2w2(t) , 



(48) 



where W 2 (t) is allowed to depend on e (although we do 
not indicate it explicitly) . Substituting into Eq. IJT7J) and 
requiring that all terms of order lower than e 4 vanish 
identically leads to the differential equation 



d_ 

Of 



W 2 (t) = V 2 (t) - D 2 (t) + i [W 2 (t), Hf(t)] . (49) 



Its general solution reads 



W 2 (t)= / duU Ht (t,u)[V 2 (u) - D 2 (u;t 2 )]U H «(u,t) 



+ U H e(t,t )B 2 U H «(t ,t), 



(50) 



where B 2 is any constant self-adjoint operator. Here we 
shall set B 2 to 0. 

The propagator UH 1 {t,to) we are looking for is ob- 
tained from Eqs. (|15|l and (|37() 

U Hl (t,t ) =T l {t)T 2 {t)U H3 {tM)T ] 2 {t a )Tl{t () ). (51) 

Substituting Eq. l(3"9"j) and neglecting Uji 3 (t,to;t 2 ) since 
it is associated with a Hamiltonian of order e 4 , we find 



(52) 



U Hl (t,t Q ) = U { H ) 1 (t,t ) + O(e 4 ), 

with the two-iteration KAM expansion 

Crg (t, t ) - T^T^UH^tM^lit^Tli^). (53) 

This expansion features two additional arbitrary times 
(t 2 and t' 2 ) and two additional arbitrary operators 
(D 2 (t 2 ;t 2 ) and B 2 ) that can be chosen so as to improve 
its accuracy. The structure of the equations and of the 
operators involved is exactly the same for the second 
KAM iteration as for the first one, and will be the same 
for any iteration. In particular, the effective perturbation 
determined at each step is always of the form of Eq. i|36|) 
since each step is just a renaming of the previous one. 
This is in contrast to the order by order methods where 
the determination of that operator requires each time 
more algebra. 

3. Summary of the time- dependent KAM algorithm 

The propagator of the perturbed Hamiltonian H\ (t) = 
Hq (t) + eVi (t) is approximated after n iterations accord- 
ing to 

U Hl (t,t ) = U$(t,t ) + O(e 2 "), (54) 



with the unitary n-iteration KAM expansion 

U$ (t, t ) = Ti (t) . . . T n (t)U H . (t, t ) 
xTt(to)...T|(t ). 



(55) 



We recall that the superexponential character of the cor- 
rection terms stems from the fact that each KAM trans- 
formation reduces the perturbation to a new effective per- 
turbation whose order is squared. One defines 



U H -(t,t ) = U Ho (t,t ) exp[-i(t - t )eDi(t ;ti)] 



x exp 



-i(t-t )e 2 D n (t ;t n ) 



(56) 



with the possibility to choose one of the following for any 
k > 1 

D fc (t;t fc )s0, (57a) 
D k {t-t k ) = t^»_ 1 (t,t fc )Vfc(tfc)l7fl S _ 1 (tfc,t). (57b) 



Each KAM transformation reads 



T k (t) = exp 



2*:- 



W k (t) 



(58) 



where 
W k (t) 



duU Hl i (t, u) [V k {u) - D k (u; t k )} U Hl i (u, t). 

(59) 

One has J7ffg(t, to) = UH (t,t ). To continue the algo- 
rithm we determine the new perturbation e 2 V n +i(t) in 
terms of the operators W n (t), V n (t) and D n (t;t n ) ob- 
tained at the preceding step 



V n+ i = 



k=l 



i k (e 2 ) fc 



k-l 



■ad k (W n ,kV n + D n ), (60) 



with the usual definition of ad(A, B) recalled in Eq. I|35|) 
and where the infinite series can be truncated to a pre- 
scribed order. This new effective perturbation has ex- 
actly the same structure at each iteration which is useful 
for applications, particularly when high orders computa- 
tions are needed. 

At each KAM iteration two arbitrary times t k and t' k 
are introduced through Eqs. I|57|l and (|59|l . As a conse- 
quence, one has the following dependence on these free 
parameters: 

UH°(t, to) = UH»(t,to;ti, ■ ■ ■ ,t n ), 

T n (t) — T n (t; ti, • • • , t n , f 1; • • • , t n ), 
W n (t) = W n (t;t ir -' ,t n ,t' ir -- ,t' n ), 
V n (t) = V n (t;t lr -- ■.. ,t' n _ 1 ). (61) 

These quantities, together with the choice of Eq. I|57a() 
or (|57b|) for the arbitrary operator D n {t n \t n ), may sig- 
nificantly affect the accuracy of the n-iteration KAM ex- 
pansion. 
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4- Comparison with the Magnus expansion 

The Magnus and KAM expansions differ in several re- 
spects. First it is remarkable that the KAM algorithm 
can be implemented in the original representation. In 
Appendix [U] we show that the result obtained for the 
KAM expansion in the interaction representation is iden- 
tical at any level of approximation. 

Most importantly of course is the superexponential 
character of the KAM expansion which manifests itself 
as of the second iteration. However, the first iteration 
of these algorithms are generally different, owing to the 
non-commutativity of the operators involved. Indeed, 
Eq. I|12|) gives for the first order Magnus expansion 

U ( H 1) i (t,t )=e-^ M ^U Ho (t,t ), (62) 

while the one-iteration KAM expansion of Eq. I|31[) reads 

Un\t to) — e ^ ie { M i(* ; *o)~ M i( t 'i ; *o)+(*i~ t )- D i( t ' t i)} 

x e -ie(t-to)Di{t;t 1 ) e -ie{M 1 (t' 1 ;t ) + (t -t' 1 )D 1 (t;t 1 )} 

xU Ho (t,t ). (63) 
We recall from Eq. lfT3)l that 

Mi(t;t ) = f duU Ho (t,u)Vx(u)U Ho (u,t). (64) 

J to 

To compare these expressions we shall cast the product 
of the exponentials of Eq. (|63|l into a single exponential 
using the Campbell-Baker-Hausdorff formula H?} 

e A e B = e c , (65a) 

where 

C = A + B + ±[A,B] + ±[A-B,[A,B]]+-.. . (65b) 

We note that the exponents of Eq. I|63|l precisely sum up 
to that of Eq. i.e., eM x (t). Hence, by Eq. l|65L) . 

these expansions differ by terms of order e 2 . In other 
words, these expansions differ through terms whose order 
is that of their remainder, which therefore enables one to 
recover precisely the same expansion up to a given order. 
To compute explicitly these terms of order e 2 and show 
that they generally do not vanish we apply the Campbell- 
Baker-Hausdorff formula twice to reduce the three expo- 
nentials of Eq. to a single one 

U^\t,t ) = e -^ M ^ t o)~4- K ^ t - t o^A)+o(<fl) 

x LWMo), (66) 

where 

*fi(t;t ,ti,ti) = {t - t' 1 )[M 1 {t-t ),D 1 {t-t 1 )] 
+ (to-i) [Mi(ti;t ),X>i(i;ti)] 
+ [M 1 {t;t ),M 1 (t' 1 ;t )}. (67) 



We recall that, according to Eq. l(2"3*|) . we choose either 
Di(t;h) = or £>i(*;*i) = t/ff (*,ii)Vi(*i)f/ ffo (ti,t). 
As a consequence Ki(t; to, t±, t[) is generally non zero 
and therefore, at this first level of approximation, the 
KAM and Magnus techniques differ by terms of order 
e 2 . However, if one chooses Di(t;ti) = together with 
t[ = to, then the one- iteration KAM expansion and the 
first order Magnus expansion coincide (this will no longer 
be true for the next levels of approximation). 

Note that for the KAM algorithm, as we discuss in the 
following section, the choice of D%(t; t\) with an arbitrary 
t\ allows one to enhance the convergence precisely by 
acting on these higher order terms (we emphasize that 
these terms appear here as higher order ones because of 
the use of the Campbell-Baker-Hausdorff formula; but at 
the level of Eq. (|63|l . each exponent is indeed of order e). 

For higher orders and iterations however the KAM al- 
gorithm is a priori expected to perform far better owing 
to both the superexponential character and the possibil- 
ity to enhance the accuracy. 



III. IMPROVING THE ACCURACY 

The formulation of the time-dependent KAM tech- 
nique presented in Sec. Ill Dl reveals the existence of sev- 
eral degrees of freedom which are at our disposal to re- 
duce the error without prior knowledge of the exact solu- 
tion: i) the choice of Eq. (|57Jl for the operators Dk(t;tk), 
and ii) the free parameters tk and t' k of Eqs. (|57bl) and 
(|59|1 . There is also a third way discussed below: iii) the 
possibility to consider another identification of the per- 
turbation and unperturbed Hamiltonian. 

The items i) and ii) also apply to the Poincare-Von 
Zeipel and the Van Vleck techniques, although to a 
smaller extent, as we describe below. 



A. Choice of Dk(t;tk) and correspondence between 
resonances and secular terms. 

Each iteration of the time-dependent Poincare-Von 
Zeipel, the time-dependent Van Vleck and the time- 
dependent KAM algorithms features an arbitrary oper- 
ator Dk(tk',tk). In the preceding section, in addition to 
the simplest case Dk(tk',tk) — 0, we suggested the choice 
Dk{tk',tk) = Vfc(ifc) where tk is an arbitrary time. 

The first iteration involves the operator D\(t; t\) which 
satisfies the same equation in the three algorithms, 
namely Eq. (|22l) . This equation is actually the general 
solution to the differential equation 

[H (t),D 1 (t;t 1 )}=i^-D 1 {t;t 1 ). (68) 
at 

This latter equation, together with Eq. I|28|) for W\(f), are 
the time-dependent generalization of the so-called coho- 
mology equations [29j considered in the stationary case. 
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The time-independent problem, i.e. the problem of 
finding a transformation T± that enables one to sim- 
plify the time-independent Hamiltonian Hi according to 
T^HiTi = Hq + eD\ + e 2 V2, is recovered when one con- 
veniently chooses T\ as time- independent. This transfor- 
mation is sometimes called contact transformation 
or level-shift transformation [28|. 

In this case all the operators, and in particular D\ and 
Wi, are time- independent and the standard cohomology 
equations are recovered 

[if 0) -Di] =0, (69a) 
Vi-Dx+i \W U H ] = 0. (69b) 

Their solutions can be determined using the following 
key property [2j|: W\ exists if and only if Hh {Di — 
Vj.) = 0, where Hh is the projector in the kernel of the 
application A i— > [A, Hq] (for an operator A acting on the 
same Hilbert space as Hq). The projector Hh applied 
on an operator A captures thus all the part B of A which 
commutes with Hq: [B, Hq] = 0. The unique solution D\ 
allowing W\ to exist and satisfying Eq. (|68|l is thus 

D 1 =Il Ho V 1 = lim I / e- im »V x e im °. (70) 

The resonances are associated with terms of V\ which 
commute with Hq. Application of Eq. (|70|l can be inter- 
preted as an averaging of V\ with respect to Hq which 
allows one to extract resonances. 

For the time- dependent problem, the general solution 
to Eq. (|28|l is given by Eq. 1|29[1 . Defining the average 

1 f* 

n_V x = lim - / dsC/ ffo (t,s)Fi(s)Z7H (s,t), (71) 

r^oo T J t _ T 

one can show the following property: if W\{t) is bounded 
for negative infinite times, then Il_(Vi — Di) = 0. This 
is satisfied by Di — Tl_Vi, the only solution compati- 
ble with Eq. (|68|l and the projector IT_. This means 
that the averaging D\ = H-V\ allows to remove secular 
terms at negative infinite times. We remark that this 
definition of the average, Eq. (|71J) . can in fact be recov- 
ered from the formal calculation of the average Hk Vi of 
Eq. I|70(l with respect to Kq = — i + Hq in an extended 
space, which includes time as a coordinate This 
gives the precise correspondence between the resonances 
of stationary problems and the secular terms of time- 
independent problems. 

In Ref. it was shown for perturbations that are 
localized in time, in a finite but possibly large interval 
ti <t < tf, that Eq. (|71|l reduces to 

n_Vi = UH (t,U)Vi(ti)U Ho (ti,t). (72) 

This is a particular solution to Eq. I|68|) corresponding 
to the choice t\ = ti in Eqs. and (|23b|l . An alternate 
definition of the average 

1 f t+T 

n+^EElim-/ d«l7fl- (t,a)Vi(s)^ ff0 (« s t), (73) 



gives a different result 

n+^i = UHofrtjKttriUH&fS), (74) 

and allows one to remove secular terms at positive infinite 
times. Generally one cannot remove simultaneously the 
secular terms at negative and positive large times. This 
shows a conceptual difference between stationary reso- 
nances and secular terms associated with perturbations 
localized in time. Furthermore, it suggests that com- 
bining both definitions in a non-trivial way gives a new 
secular term that could improve the convergence of the 
algorithm. This is achieved by the general solution to Eq. 
iJHSJ, Eq. H22fl . written with the perturbation evaluated 
at a free time t x as the arbitrary operator [cf. Eq. 123b(l ] 

Di(t;ti) ee C/if (Mi)Fi(ti)[/tf (ti,i). (75) 

The free t\ can then be chosen as we describe below to 
minimize the remainder obtained at the first iteration of 
the perturbative algorithm. 

The next iterations also offer the possibil- 
ity to choose for Dk(t;tk), in particular or 
UHi^fatkMit^Uffe ^t) similarly to Eq. {Z3- 
Note however that for the Poincare-Von Zeipel and the 
Van Vleck techniques, the choice corresponding to Vfe(tfe) 
can only be made once (for an arbitrary value of k). As 
explained in Appendixes A and B, respectively, this is 
due to the order by order character of these techniques. 

B. Enhancing the convergence 

After one iteration of the KAM algorithm one deduces 
from Eqs. I|15f) and (|17|l an exact expression for the full 
propagator 

Uh x (t, t ) — Ti(t; ti,t[)UHi(t, t'{; ti)Un 2 (t, t ; 

x UHvit'lMijTtita-MA), (76) 

where t\, t[ and t'{ are arbitrary times [cf. Eq. I)32fl ]. The 
propagator Uu 2 (t, to]ti,t'i,t1) of the Hamiltonian given 
in Eq. (|18|l is associated with a second order generator 

U R2 {tMhAJl) = e-^Wo.'^.'"). (77) 

To obtain the one-iteration KAM expansion Ug\t,to) 
we neglected this propagator, replacing it by the iden- 
tity in the above product. Obviously, the closer 
Un 2 (t, to; ti, t[, t'l) is to the identity, the smaller the cor- 
rection terms are, i.e., the more accurate the approxima- 
tion is. We can improve this accuracy if we can make that 
propagator closer to the identity, or equivalcntly, its gen- 
erator closer to zero. The distance is defined through the 
norm \ \A\ \ = max||^|| =1 | \Axj)\ \ with ip in the Hilbert space 
of the problem. For an Hermitian matrix, this norm is 
the largest absolute value of its eigenvalues. 
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We calculate this generator e 2 G2(£; to, ti, t' 1; t'i) by 
solving the Schrodinger equation with the Hamilto- 
nian of Eq. 1)18(1 in the form of an exponential using 
Eq. I|tj6|) . This is a time-dependent problem with a 
zero unperturbed Hamiltonian and whose perturbation is 
e 2 UH-i(t'{,t;t 1 )V 2 (t;t 1 ,t' 1 )UHi(t,t'{;t 1 ). Hence, we evalu- 
ate the lowest order contribution to G*2(£; to 7 *i , t[, t") as 

G^ 2) (t; t , ti, t'i) = / du L%e(io,-u;ti)V2(u;ii,ii) 

x U H1 (u,t ;h), (78) 

where we set t'[ = t (its precise value is not relevant 

since the one-iteration KAM expansion U^(t,to;tx,t'x) 
is independent of the parameter t'{). It is this operator 

G 2 (t; to, tx, t'i) that has to remain small for the algo- 
rithm to converge |35j . Having the arbitrary times tx 
and t'i at our disposal, we can actually enhance the con- 
vergence of the algorithm by minimizing the norm of this 
operator with respect to these free parameters. 

Similarly, the n-iteration KAM expansion of Eq. I|55|) 
can be optimized by minimizing the norm of the follow- 
ing operator with respect to one or several of the free 
parameters f 1, • • • , t n , i' l7 • • • ,t' n 

G%g\t)= / duU H *(t ,u)V n+ i(u)U H e(u,t ) 

J t 

= @n+X " " j *ru *1> — , O- (79) 

The dependence of t ) and V^+i(t) on these pa- 

rameters is given in Eq. 

It turns out, as will be illustrated in Sec. IIV Bl that 
modifying the parameters tk and/or t' k can improve the 
accuracy by more than one order of magnitude already 
for k = 1. 



C. KAM expansion with another identification of 
the unperturbed Hamiltonian 

The perturbed Hamiltonian can be written as 
Hx (t) = H (t) +eDx(t;tx) + eVx(t)-eDx{t;tx) 

V v ' " * ' 

= H' (t;tx) + eV{(t;tx), (80) 

where Dx(t;tx) = Uh (t, t x)Vx ( tx)UH (tx , t) with tx ar- 
bitrary [cf. Eqs. li-'L'l) and (|23b(l ] . The propagator asso- 
ciated with H'o(t]tx) can always be determined since by 
Eqs. (|T§j) and lf251) one has 

H' Q {t;tx) = Hl(t;tx), (81a) 
U Ho (t,t ;tx) = U H! (t,to;tx). (81b) 

We now apply the KAM algorithm exactly as summa- 
rized in Sec. Ill D3l i.e. with the same definitions for all 
the operators involved in the expansion, but with the 



identification of Eq. (jSDJ. This decomposition has the 
property that V{(tx;tx) = which implies by Eq. (|23b|) 
that Di(t; tx) — 0. The free parameter tx is therefore in- 
troduced here through Eq. l|81bH . One arrives at a KAM 
expansion which is still of the form of Eq. (|55|l but may 
significantly differ from that resulting from the conven- 
tional decomposition. 

We emphasize that the possibility to consider the iden- 
tification of Eq. (|SU|) as a new starting point for a pertur- 
bative treatment is specific to the KAM technique which 
is not an order by order method, contrary to the Magnus, 
the Poincare-Von Zeipel and the Van Vleck algorithms. 

IV. BEYOND THE SUDDEN APPROXIMATION 
A. Preliminaries 

We consider a system described by the Hamiltonian 
H (autonomous or not). It is perturbed by a time- 
dependent Hamiltonian V(t) whose characteristic dura- 
tion is r. This latter quantity is the time during which 
the interaction differs significantly from zero, and not the 
full duration of the interaction, which may be large but 
finite. Here we define r as twice the full width at half 
maximum, having in mind a perturbation which presents 
a time-dependent envelope. We assume that the pertur- 
bation V(t) satisfies 

[V(t),V(t )]=0 Vt,t , (82) 

which is realized in many situations of physical interest. 
The propagator of the perturbed system evolves accord- 
ing to the Schrodinger equation 

z^U(t,t ) = {H+V(t)}U(t,t ), (83) 

with U(to,to) = 1, the identity operator on the appropri- 
ate Hilbert space EL We define a dimensionless time t and 
dimensionless operators H 7 V(t) and U(t, t ) through 

t = rt, 
H = HujH, 

V(t) = -V(t), 

T 

U(t,t ) = [/(Mo), (84) 

where u> is some characteristic frequency of H . In dimen- 
sionless units Eq. H83fl becomes 

i?rU{t 1 t ) = {V{t) + eH}U{t 1 t ), (85) 
at 

where we define a sudden parameter e = lot . 

The sudden or impulsive regime corresponds to the 
limit e — > 0. Our aim is to obtain a perturbative expan- 
sion for the evolution operator U(t, to) of the perturbed 
Hamiltonian V(t) + eH beyond the sudden regime. To 
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this aim we identify the original perturbation V(t) as the 
unperturbed Hamiltonian H (t) and the original unper- 
turbed Hamiltonian H as the perturbation V\ 



H (t) = V(t), 
Vi = H. 



(86a) 
(86b) 



Note that we need only consider a finite interval of time 
as V(t) is localized in time. By virtue of Eq. I|82|l. the 
propagator for e = reads 



Uh (t, t ) = exp < -i / V(u) du 



(87) 



We shall follow this approach below and consider the 
various pcrturbative schemes described in Sec. ^ in the 
case of two-level systems driven by short pulses. 



B. Illustration on pulse-driven two- level systems 

Our purpose is to compare the various algorithms and 
investigate the convergence enhancement as well as to 
show that the unitary time-dependent KAM theory is 
well suited to study regimes beyond the impulsive or sud- 
den limit. 

In the notations of the preliminaries, we take H = fnoa^ 
and V(t) = -Cl(t)ai where fi(t) is a pulse shape function 
and a k are the Pauli matrices 



1 

1 



cr 2 



-i 

1 



CT3 



1 

-1 



Hence, the operators defined in Eq. (|5BT) are here H (t) 
Q(t)cri and V\ = 03. The Schrodinger equation reads 

d 

i-^U Hl (t ) t ) = [n(t)(7i + aT 3 ]l7ff 1 (t,*o), 
with Uh-i (to, to) = 1c 2 - For e = its solution is 



UH (t,t ) = e 



-iA(t)<Ji 



(89) 



where A(t) = J t Sl(u) du. The pulse area A = A(oo) is a 
dimensionless parameter that can be fixed independently 
of the sudden parameter e (= lot) that we take here as 
the perturbative parameter. This allows us, in particu- 
lar, to treat large nonperturbative areas for short pulse 
durations. 

The Magnus expansion (t, t ) is given by Eqs. i|12|) - 
(THall . 

The KAM expansion is obtained from Eqs. l(K5 |l -l|60 jl . 
We distinguish three types of KAM expansions, reflecting 
the choices discussed in Sees. Ill Dl IIII Al and IIII Cl 

Type A: each iteration k involves the operator 
D k (t k ;t k ) =0 [cf. Eq. J5H)]. 

Type B: D k (t k ;t k ) = V k (t k ) for all k [cf. Eq. J57bj]. 



Type C: the unperturbed Hamiltonian is defined as 
2*5 (t; *0 = H (t) + eUHofatJViitjUHoitut) [cf. 
Eq. (JHOJl]. 

For the type A n-iteration expansion one has n free times 



^Hi (t,to) — U^(t,to] t' x , ■ ■ ■ ,t' n ), 
while for the types B and C one has 2n such parameters 

^Hi (*>*o) — (t,to;tx, ■ ■ ■ ,t n ,t[, ■ ■ ■ ,t' n ). 

In the case of two-level systems, the infinite series of 
Eq. (|rJU|) for the new effective KAM perturbation V n+ i(t) 
can be cast into the form [1J| 



V n+ i = [W n , a n V n + b n D n ] 

+ e n [W n , [W n , CnV n + d n D n }} 



(90) 



where 



f 2 A 2 

71 71 



e 2 X 2 

n y n 



_ (n^n cos e„A n — sin e n X n , _ ., i<\-\\ 



with e„ = e 2 " 1 and \ n (t) = \J — det W n (t). 

For given e and A, the error A„ at the end of the pulse 
between the numerical solution of the Schrodinger equa- 
tion and the result obtained after n iterations is defined 
as 



/\n^\\U Hl (tf,U) - U { ^(t f ,U)\[ 



(92) 



We also define the error 5 n in the transition probability 
from the lower state |— ) to the upper state |+) 



<+l^(*/,*i)|-> -\(+\U Hl {t h U)\-)[\ (93) 



We consider the following dimensionless pulse shape 
between the dimensionless time U — and tf = 1 



„, s ( 2Asin 2 (irt) < t < 1, 
^ ' I elsewhere. 



(94) 



Figure 1 displays the common logarithm of the error 
Ai and the error 5\ as a function of the pulse area A 
for the Magnus expansion and the three types of one- 
iteration KAM expansions in the case t\ = t' y = and 
e = 0.5. The errors Ai and Si globally decrease when A 
increases. This is expected on the basis of Eqs. (I88|l and 
(19411 as the relative importance of the perturbation then 
decreases. One also observes marked oscillations in Ai 
and 5i with a pseudo-period of n. This stems from the 
form of the unperturbed propagator [cf. Eq. and 
the fact that it always appears twice, in particular in the 
operator Gv, (t;to) of Eq. (|78|l which controls the error 
after one iteration. 
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Let us recall that the one-iteration KAM expansion of 
type A coincides with the first order Magnus expansion 
for t'i = 0. It is seen, by both measures of the error, 
that each of the KAM expansions can perform better 
than the other ones on some intervals of A. Hence, in 
order to establish a fair comparison, we shall consider 
the particular value A = 1 where the first order Magnus 
expansion and the one-iteration KAM expansion of type 
B yield essentially the same error Ai for e = 0.5. This 
remains true for all values of e up to 2 as can be seen from 
Fig. 2 which depicts the errors Ai and Si as a function 
of e for A = 1. 

In Fig. 2 we also present the one-iteration KAM expan- 
sion of type B that is optimized by choosing t\ — 0.5, t\ = 
0.22. We see that the error Ai is reduced (with respect to 
the comparable Magnus and non-optimized type B KAM 
expansions) by more than one order of magnitude up to 
values of e equal to unity. The error 8\ on the transition 
probability is also considerably reduced. Notice that the 
values of for the Magnus and non-optimized type B 
KAM expansions differ while the values of Ai are indis- 
tinguishable. For comparison, we also consider the (non- 
unitary) Dyson expansion [23 ■ Recall that it is obtained 
by repeated use of the integral form of the Schrodinger 
equation in the interaction representation with respect to 
H {t) 



U H i(t,t ;to) = l-i / duHl^t^UHiiu^to) 



= l- i duH[(u;t ) 



1 - i I dvH{{v;t ) 

to 



(95) 



where H[(t;to), given by Eq. l(2*c|) . contains a prefactor 
e. One then returns to the original representation with 
Eq. (|2a|l . For the value of A considered in Fig. 2, the 
Dyson expansion yields the largest error Ai whereas its 
error on the transition probability is rather small. 

In Fig. 3 we plot g 2 defined as the largest abso- 
lute value of the eigenvalues of the Hermitian matrix 
*l>*'l) given by Eq. (JJSJ. This quantity which 
controls the error after one iteration is represented as a 
function of t\ and t[ for the KAM expansion of type B. 
By minimizing g 2l which is the norm of this operator, 
with respect to the free parameters t\ and t' 1; one re- 
duces the error without having to determine the exact 
solution. Notice that (t\ — 0, t[ — 0) is a local maxi- 
mum of g 2 whereas (t\ — 0.5, t[ — 0.5) is a saddle point. 
The point (t\ = 0.5,^ » 0.22) corresponds to a mini- 
mum. The symmetry of Fig. 3 results from the pulse of 
Eq. 194|) being symmetric. 

Figure 4 displays the error Ai and the eigenvalue g 2 as 
a function of t' x for the three types of one-iteration KAM 
expansions in the case A = 1, e = 0.5. The value of ti 
is chosen so that the optimum can be reached: t\ = 0.5 
for the type B and t\ = 0.7 for the type C (recall that 



the type A features no t\). One sees that the error can 
be reduced by more than one order of magnitude for the 
type B or C, and about half an order of magnitude for 
the type A after a single KAM iteration. It is also seen 
that the eigenvalue g 2 is a very accurate estimation of the 
error Ai, which enables one to locate the optimal values 
of the free parameters. Note that the first order Magnus 
expansion corresponds to the particular value t\ = (i.e. 
the nonoptimized case) of the one-iteration type A KAM 
expansion. 

We now turn to the next level of approximation for the 
Magnus expansion and the KAM expansion of type B. 
Recall that the first order and one-iteration expansions 
of these schemes yield comparable errors Ai for A = 
1. Figure 5 shows that the (nonoptimized) two-iteration 
KAM expansion performs better than the second order 
Magnus expansion by one to two orders of magnitude for 
A2. The error S 2 on the transition probability is also 
much smaller for the KAM expansion. 

It is worth noting from the comparisons of Figs. 2 and 
5 that the error Ai for the optimized one-iteration KAM 
expansion is comparable to the error A2 for the nonop- 
timized two-iteration KAM expansion of type B. This 
conclusion is not restricted to the type B and can be un- 
derstood on the basis of the Campbell-Baker-Hausdorff 
formula as discussed in Sec. lIID4l 

If one optimizes the type B KAM expansion by choos- 
ing t\ = 0.5, t[ = 0.22 (as determined above from Fig. 3) 
and t 2 — 0.66, t' 2 = 0.8, one gains another one to two 
orders of magnitude on A2 . 

From Fig. 5 one also deduces that the Dyson approach 
is not applicable in this context as the second order per- 
forms worse than the first order by both measures of the 
error. Notice that the transition probability predicted by 
the Dyson technique diverges, as is well-known, by lack 
of unitarity. In other words, the Dyson expansion does 
not allow one to refine the results of Fig. 2. 

The two-iteration KAM expansion involves the opera- 
tor V 2 (t) which is given by Eq. (|34|l or l|36|) as an infinite 
series of commutators. For two-level systems, this series 
can be computed explicitly and results in Eq. I|90|) with 
ri = 1. It is remarkable that the coefficient a±, bx, c\ 
and d\ are well-defined for all values of e, even larger 
than unity. In Fig. 5 we consider the cases where V 2 (t) is 
truncated to two commutators [i. e. the term k = 1 of or- 
der e° and the term k — 2 of order e in Eq. I|36ll ]. and four 
commutators (k = 1, • • • , 4). Note that this amounts to 
approximate the coefficients of Eq. I|91|l by polynomials 
(of order, respectively, 2 and 4) in eAi(t). The case of 
two commutators performs better than the second-order 
Magnus expansion which is not surprising as it contains 
all the terms of order e 3 . However, it performs worse than 
the KAM expansion with the infinite series by up to one 
order of magnitude. The case of four commutators is 
very close to the exact two-iteration KAM expansion for 
values of e up to unity. The convergence with the number 
of commutators involved is indeed very fast. It has to be 
remarked that including commutators of higher orders 
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in a well-defined manner as in the KAM algorithm is the 
main difference with the Magnus expansion (or any order 
by order expansion) and necessary to achieve superexpo- 
nentiality. From a practical point of view the inclusion of 
more than one commutator is both convenient and highly 
profitable. 

The fact that for the KAM expansion the order of the 
remainder be the exponential of an exponential (e 2 ) 
implies that the accuracy is increased in a dramatic way 
with the number n of iterations. In addition, the ac- 
curacy of a given n-iteration KAM expansion can still 
be enhanced, as illustrated here, and allows one to get 
closer to the next iteration of the algorithm by including 
through the simple variation of the parameters t n and t' n 
more appropriate higher order terms. 
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FIG. 1: Comparison of the first order Magnus expansion 
(dots), the one- iteration type B KAM expansion (solid line) 
and the one-iteration type C KAM expansion (dashed line) for 
e = 0.5 and ti = t\ = 0: (a) common logarithm of the error 
Ai and (b) error Si as a function of A. The one-iteration type 
A KAM expansion coincides in the case t'i = with the first 
order Magnus expansion. All quantities are dimensionless. 
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FIG. 2: Comparison of the first order Magnus expansion 
(dots), the first order Dyson expansion (circles), the one- 
iteration type B KAM expansion with ti = t'i = (solid 
line) and the optimized one-iteration type B KAM expansion 
with ti = 0.5, ti = 0.22 (dash-dot line) for A = 1: (a) com- 
mon logarithm of the error Ai and (b) error Si as a function 
of e. All quantities are dimensionless. 



the other methods, has been established owing to its su- 
perexponential character and the accuracy optimization. 
We have also shown that the Magnus expansion is recov- 
ered as a special case of the time-dependent Poincare-Von 
Zeipel expansion (whose time-independent version coin- 
cides with the Rayleigh-Schrodinger expansion [25|). 

The possibility to enhance the accuracy of a given level 
of approximation stems from the free parameters and the 
free operators that appear naturally in the formulation 
presented here. It allows one to significantly reduce the 
error with respect to the exact solution without its knowl- 
edge by the minimization of an eigenvalue. 

The above considerations, illustrated here on a pulse- 
driven two-level system, are straightforwardly applied to 
more involved problems (see Ref. |30| for an application 
to the orientation and alignment of molecules). 



V. CONCLUSIONS 

We have formulated perturbation theory in operator 
form |28j for time-dependent problems localized in time 
directly in the original Hilbert space by unitarily trans- 
forming the evolution operator. 

We have compared formally and numerically various 
unitary perturbative schemes. The superiority of the 
KAM technique over the Magnus expansion, as well as 
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FIG. 3: Eigenvalue for the KAM expansion of type B as 
a function of ti and t[, for A = 1 and e = 0.5. All quantities 
are dimensionless. 



APPENDIX A: TIME-DEPENDENT 
POINCARE-VON ZEIPEL EXPANSION 

The time-dependent Poincare-Von Zeipel technique 
amounts to construct a unitary operator T(t) which 
transforms the propagator UH 1 (t,to) into a propagator 
UH*(t,to) according to 

Ti(t)U Hl (t,t )T(t ) = U m (t,t ), (Al) 

where UH«(t,to) is associated with an effective Hamilto- 
nian H e (t) containing contributions to every order in e 



H e (t) = H (t) + ^2e k D k (t;t k ), (A2) 
fe=i 

This propagator can be expressed in terms of Uh (t, io) 
and unitary operators related to the partial Hamiltonians 
e k D k (t;t k ) 

U H *(t,t ) = U Ho {t 7 t )exp [-i(t - t )e-Di(*o;*i)] 

x ■ ■ • exp [-i(t - t )e k D k (t ; t k )] • ■ • , (A3) 

provided the operator D k (t; t k ) are defined as 

D k (t;t k ) = U H ._ 1 (t,t k )D k (t k ]t h )U Ht _ 1 (t k ,t). (A4) 

Here D k (t k ;t k ) is arbitrary but strictly of order e° (in 
contrast to the KAM case) and we set 

U H °(t,t ) = U Ho (t,t )exp [-i(t - t Q )eD 1 (t ;ti)} 

x • ■ • exp [-i(t - t )e n D n (t ; t n )} . (A5) 
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l 1 

FIG. 4: Comparison of the one-iteration type A KAM ex- 
pansion (dots), the one- iteration type B KAM expansion with 
ti = 0.5 (solid line) and the one-iteration type C KAM ex- 
pansion with ti — 0.7 (dash-dot line) for A = 1 and e = 0.5: 
(a) common logarithm of the error Ai and (b) eigenvalue g2 
as a function of t[. All quantities are dimensionless. 

The generator of the Poincare-Von Zeipel transformation 
T(t) is written as a power series of e-independent opera- 
tors W k (t) 

T(t) ee exp ^-ijt^W^ . (A6) 

These operators satisfy the differential equations 

^-W k (t) = V k (t) - D k (U t k ) + i [W k (t), H (t)} , (A7) 
at 

where the expression for V k (t) cannot be given in a sim- 
ple form for arbitrary k and requires increasing algebra. 
The situation is analogous for the Magnus expansion or 
the Van Vleck expansion given below, and in contrast to 
the KAM expansion where the new effective perturbation 
has exactly the same form at each step [cf. Eq. (|60(l ]. 
The general solution to Eq. I|A7J| reads (up to a term 
UH {t,to)B k UH (to,t) with B k any constant self-adjoint 
operator) 

W k {t)= f duU Ho (t,u)[V k (u)-D k {u)}U Ha (u,t), (A8) 

where t' k is arbitrary. Notice that this expression involves 
the unperturbed propagator unlike the KAM analogue 
which features an effective propagator [cf. Eq. (|59|) ]. 
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FIG. 5: Comparison of the second order Magnus expansion 
(dots), the second order Dyson expansion (circles), the two- 
iteration type B KAM expansion it\ — t' x = t 2 = t' 2 — 0) 
with the infinite series of commutators (solid line), two com- 
mutators (dashed line), four commutators (squares) and the 
optimized (ti = 0.5, t'i = 0.22, t 2 = 0.66, t 2 = 0.8) two- 
iteration type B KAM expansion (dash-dot line) for A = 1: 
(a) common logarithm of the error A2 and (b) error Si as a 
function of e. All quantities are dimensionless. 



For the operator Dk(t;tk) of Eq. (|A4|I to be strictly 
of order e°, we are no longer entitled to choose 
U H i_ 1 (t, tk)V k {t k )U H -- k _ 1 (tk,t) for any k if that choice has 
already been made for a lower value of k. This choice can 
only be made once (for an arbitrary value of k denoted 
v) and implies to take Dk(t;tk) = for the other values 
of k. 

We then construct the n-th order unitary approxima- 
tion to T(t) 



fc=i 



(A9) 



and by Eqs. IjAlfl and i|A3(l obtain the n-th order 
Poincare-Von Zeipel approximation to the propagator 
U Hl (t,t ) 

U$(t,t )=T n (t)U Ho (t,t ) 

x exp [-i(t - t Q )e v D v (t ; t v )] T*(to), (A10) 

where v is any integer between 1 and n for which 
we have the possibility to choose either D v (t ;t v ) = 
UH (to,t v )V v {t v )U H „(tv,to) or D v (t ;t v ) =0. 



It is interesting to note that Eq. I|A10|) precisely re- 
duces to the Magnus expansion if we take D v (t;t v ) = 
and t' k — tg for all k. In the general case there are up to 
n + 1 free parameters 



T n (t) 
W n (t) 



W n (t;t v ,t' V ' 



t' ) 



(AH) 



APPENDIX B: TIME-DEPENDENT VAN VLECK 
METHOD 

The time-dependent Van Vleck technique is an order 
by order method where instead of transforming the origi- 
nal propagator into a final new propagator in a single step 
as in the Poincare-Von Zeipel algorithm, one achieves this 
goal iteratively through a series of transformations T k (t) 
which reduce the size of correction terms from e k to e k+1 . 
This results in the n-th order Van Vleck expansion 

(t, t ) = Ti(t) ■ • ■ T n (t)U Ho (t, t ) 

x exp [-i(t - t )e v D v (t ;t v )] T+(t ) • ••T-ftto), (Bl) 

where D v (t ;t v ) = U Ho (to,t v )V v (t v )U Ho (t v ,to) or 
D v (to] t v ) =0, with v an integer between 1 and n. Here 
we define 

-ie k W k (t) 



T k {t) 



(B2) 



with Wk(t) constructed as in Eq. (JA8J). There are up 
to n + 1 free parameters entering the n-th order time- 
dependent Van Vleck expansion. 

The Van Vleck and Poincare-Von Zeipel techniques dif- 
fer because of the product of exponentials of Wk(t) ap- 
pearing in Eq. (|F31|) instead of the single exponential of 
a sum of Wk(t) in Eq. IjAlOp . As discussed above on 
the basis of the Campbell-Bakcr-Hausdorff formula this 
means that these algorithms differ at orders higher than 
the prescribed order. 



APPENDIX C: KAM EXPANSION IN THE 
INTERACTION REPRESENTATION 

In the framework of a perturbation theory, we stressed 
111 Sec. IITaI that the Ma gnus expansion had to be derived 
in the interaction representation. In the section above, 
the KAM algorithm was applied in the original repre- 
sentation. We show here that going to the interaction 
representation, applying the KAM technique and com- 
ing back to the original representation yields identically 
the same expansions as above despite the truncation at 
any finite order. 

Given an Hamiltonian Hi (t) = Ho (t) + eV\ (t) and the 
propagator {/# (i, to) we consider the interaction repre- 
sentation with respect to Ho(t). From Eq. (J2J one de- 
duces that H[(t; s) = H (t; s) + eVftt; s) with 

H (t;s)= 0, (Cla) 
Vl(t;s) =U Ho (s,t)Vi{t)U Ho {t,s). (Clb) 
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Furthermore, the unperturbed propagator is trivial 

U Hh (t,t ;s) = t. (Clc) 

It is worth pointing out that Eq. (|Clb|) suggests that 
eVi(t; s) is considered small with respect to H^t; s) = 0. 
As a matter of fact, by virtue of Eq. I|2bjl . it is with 
respect to i? (t; s) — i that eVl(t; s) is considered small 
(in a technical sense we need not specify here). 

Applying the KAM algorithm with the identifications 
of Eqs. (|Clbj) and (|ClcJl leads to the expansion 

U$ (t, t ; a) = Ti (t; s) ■ ■ ■ T„ (t; s) U H « (t , t ; s) 

xT^t ;s)...T}(t Q ;s). (C2) 



Returning to the original representation with the help 
of Eq. (|2a|l one obtains a perturbative expansion for 
UH 1 (t,to) that coincides exactly with the expansion 

U-y^ (i, to) obtained directly in this representation 

U Ho (t,s)U { ^(t,t ;s)U Ho (s,to) = U$(t,to). (C3) 

This stems from the following identities readily derived 
on the basis of Eq. (|10|l and valid for any k 

U Ho {t,s)T k {t;s)U„ (s,t) = r fc (t), (C4a) 
U Ho (t,s)U Hf (t,t Q ;s)U Ho (s,t ) = U H »(t,t ). (C4b) 
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